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Abstract 

Assuming that the spin a of the black hole, presumably located 
at the core of the active galactic nuclei Messier 87 (M87), takes the 
value which maximises the ergospheric volume of the Kerr spacetime, 
we find results compatible with the recent observations obtained by 
high resolution interferometry on the origin of the jet, which would be 
located inside the innermost stable circular orbit diameter. Moreover, 
we find a flow of unbound geodesies issued from the ergoregion able 
to frame the best fits at large scale recently obtained for describing 
the observed profile of the relativistic jet launched from this central 
engine. 



1 The jet model 

The arrival of high resolution observations, of the order of microseconds of 
arc, will allow us soon to observe the origin of extragalactic jets and specially 
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the less distant ones like the M87. These data should tip the balance in favor 
of one of the two main kinds os models. The first ones situate the origin of 
the jet close to the interior of the accretion disk, near the innermost stable 
circular orbit (ISCO) (Blandford & Znajek 1977). While the second ones 
consider that the jet originates inside the ergosphere of the Kerr black hole 
(BH) situated at the centre of the galaxy, for instance proposed in (Williams 
1995, 2004) and (Gariel et al. 2010; Pacheco et al. 2012). 

The model (Gariel et al. 2010; Pacheco et al. 2012) is based on the 
unbound particle geodesies leaving the Kerr BH ergosphere. These particles 
are produced by a Penrose process (Penrose 1969) inside the ergosphere and, 
while ejected, follow asymptotic geodesies parallel to the axis of symmetry z 
at a radial distance 

r q i 1/2 

P = Pl = a [ 1 + a^-l) \ ' (1) 

where Q is the Carter constant, a the BH spin and E the particle energy. 
There is a limited number of discrete values p\ for which the energy E of 
the particles tends to infinity, while near the neighborhood of each of these 
values the energy diminishes steeply. This permits modeling a thin jet highly 
energetic and collimated. 



2 Maximal ergospheric volume 



Admitting that the Penrose process is the source of production of particles 
building the jet, it is reasonable to suppose that its frequency increases pro- 
portionally to the ergospheric volume, which is the only region where it can 
take place. Hence, the larger is the ergospheric volume the larger is the 
amount of particles produced to take part of the jet, and the more powerful 
becomes the jet. 

The trace of the ergosphere in the plane (p, z) is a closed curve, with 
equation 



(2) 



depending only on the parameter a, the Kerr BH angular momentum per 
unit mass and where we normalized its mass M — 1. The surface S inside 
this curve can be calculated from 
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which, after integration, produces 



S= — r 
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We can plot the function S(a), figure 1, for a G [0, 1] and obtain the value of 
a where S is maximal: 

a = 0.877004. (5) 

The corresponding trace of the ergosphere, (jSJ), for a given by (jSJ), generating 
the maximal ergospheric volume, is plotted in figure 2. 




Figure 1: Plot of the ergoregion projection S in the (p, z) plane as a function of 
BH spin a (see ©). Its maximum S max = 0.421385 is for a = 0.877004. 



3 Jet produced by M87 

The first observation of an extragalactic jet dates back to the year 1918 
and concerned the one produced by M87, being very long and powerful. 
Currently, the M87 nucleus distance is estimated at 16.7±0.6 Mpc (Blakeslee 
et al. 2009) and the mass of its central BH to be M = (6.2 ± 0.4) x 10 9 M sun 
(Gebhardt et al. 2011). Its relativistic jet extends for hundreds of kpcs before 
ending in large radiolobes. 
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Figure 2: Plot of the ergosphere trace, in the (p, z) plane, limiting the maximal 
ergoregion obtained in figure 1 for a = 0.877004. We also plot the beginning of 
some unbound geodesies with high energy, E = 10 6 , from this maximal ergosphere, 
of interest for modelling the M87 jet, as we can see in sections 4 and 5. The 
geodesies asymptotic to p\ = 3.1637 is obtained for the i.e. (pi = 0.565256, Z{ = 
—0.25) (see section 4), and a set of imperfectly collimated near geodesies framed by 
two geodesies, and "external" one with i.e. (pi = 0.7, Z{ = —0.4), and an "internal" 
one with i.e. (pi = 0.85, z% = —0.15) (see section 5). 

In the framework of the model discussed in (Gariel et al. 2010; Pacheco 
et al. 2012), it is natural, according to the results of the section 2, to expect 
that the most energetic jets produced by a Kerr BH of mass M are obtained 
for a given by (jSJ). It is natural to assume too that the jet produced by M87 
fits into this model. We call the attention that this value for a is consistent 
with the few evaluations of the M87 BH spin obtained to date, like a > 0.65 
in (Wang et al. 2008), and more recently a > 0.8 in (Li et al. 2009). 

Hence we look for particles moving along unbounded geodesies with high 
energy, we take E = 10 6 , produced in the ergoregion, which means, in 2D, 
produced between the ergosphere and the event horizon situated on the z 
axis. These particles move asymptotically parallel to the z axis with p — > 
pi when z — > oo. In fact, there are other geodesies, near these perfectly 
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collimated ones, that asymptotically behave p — > oo when z — > oo, however 
they are not very divergent, i.e., they satisfy p <C z for a long distance along 
the axis. 

Now we need to determine the parameter p\. In order to do this we 
consider the recently obtained high resolution data (Doeleman et al. 2012). 

4 Jet near the core of M87 

The data presented in (Doeleman et al. 2012) opens a new era of resolutions 
never before achieved, which are of the order tens of microseconds of arc. 
Hopefully soon it can be expected resolutions even better, of microseconds 
of arc. With these observations the models considering jets produced by 
the accretion disk near the ISCO (Blandford & Znajek 1977) and models 
with jets leaving the ergosphere (Williams 1995, 2004; Gariel et al. 2010; 
Pacheco et al. 2012) will be discriminated. Among these models the one 
that is based only on gravity to explain the genesis of jets considers that 
the Kerr spacetime has the essential features to explain their formation, 
acceleration and collimation (Gariel et al. 2010; Pacheco et al. 2012). Other 
fields, excepting the gravitational, like the magnetic field are no more that a 
complement for their formation. 

Through interferometry by using the VLBI array and by choosing a fre- 
quency of the order 229 Ghz (or wave length 1.3 mm) allowing a good trans- 
parency up to the core3 of the galaxy and with unprecedent resolution of the 
order of 40 ± 1.8// as corresponding, for a distance 16.7±0.6 Mpc, to a length 
of the order of 5.5 rs, Doeleman et al. (Doeleman et al. 2012) evaluated the 
maximal diameter of the emerging jet at 

r jet Aw = 5.5 ± 0.4r 5 , (6) 

where rs = 2M is the Schwarzschild BH radius. 

Then, by admitting the validity of the model by Blandford and Znajek 
(Blandford & Znajek 1977), they assume that the jet starts at the ISCO of 

1 We call the "core" the region where the jet is produced, which means extending 
maximally up to the radius of the ISCO and, consequently, including the ergosphere. This 
"core" has an altitude z of the order of M. We precise this point since it is not always 
understood in this way. For instance in (Hada et al. 2011) it is defined centred transverse to 
the axis z. Hence there are "cores" for each altitude z according to its observed frequency 
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the accretion disk and put 

rjetApp = TiscOApp- (7) 

This assumption implies too that the origin of this observation is on the 
plane of ISCO, i.e. the equatorial plane z = 0. The apparent radius of ISCO, 
plotted in (Doeleman et al. 2012) as a function of the spin a (see figure 3 
in (Doeleman et al. 2012)), is the intrinsic radius corrected by the effect 
of gravitational lensing caused by its proximity to the core endowed with 
a Kerr BH. In this figure, the horizontal line representing rj et A PP intersects 
the prograde radius of ISCO, which allows to evaluate the spin of the BH 
as a ~ 0.625 ± 0.25 and to show that only the prograde part is viable. We 
observe that the so obtained value of a is inside the limit of (Wang et al. 
2008), but outside the limit of (Li et al. 2009). 

Now we interpret the same result in the light of the model given by (Gariel 
et al. 2010; Pacheco et al. 2012). Instead of considering the "apparent" 
distances, we deduct the lensing and consider the "intrinsic" distances 

Tjetlnt ~ r je tApp ~ 1-3, (8) 

by observing that the curve riscoint(a), figure 3, roughly is translated down- 
wardly about 1.3 from the ri S coA PP {o) given in figure 3 in (Doeleman et al. 
2012). Writing these quantities in terms of Weyl coordinates, given by (1) in 
(Pacheco et al. 2012), we obtain 

Pjetlnt = [(r je tmt - l) 2 - (1 - a 2 )] 1/2 = [(r jetApp - 2.3) 2 - (1 - a 2 )} 1 ' 2 . (9) 

Now we can plot the curves Pj e tint(a) for values in the range Tj et A PP £ 
[5.1,5.9], which are given in figure 4. In particular, from the central curve, 
corresponding to the apparent radius rj et A PP = 5.5, we deduce for the maxi- 
mal volume of the ergoregion ([5]) 

Pjetlnt = 3.1637. (10) 

We abandon the hypothesis ([7]) assumed by [?] and we replace it by 

Pjetlnt = Pi- (11) 

This allows us to obtain the unbound pi-asymptotic geodesic, (see the lower 
curve, asymptotically parallel to the z-axis, in figure 5), with parameters 



6 



Pisco 




Figure 3: Plot of the Kerr prograde and retrograde ISCO for Pisco as function 
of the BH spin a obtained from the transformation formula © in the equatorial 
plane 9 = vr/2, with the standard r IS co = 3 + Z 2 ± [(3 - Zi)(3 + Z 1 + 2Z 2 )] 1/2 , 
with Zi = 1 + (1 - a 2 ) 1 / 3 [(l + a 2 ) 1 / 3 + (1 - a 2 ) 1 / 3 } and Z 2 = (3a 2 + Z 2 ) 1 / 2 . 

a = 0.87705, E = 10 6 and p\ = 3.1637 after numerical integration of (21) 
given by (Gariel et al. 2010) and initial conditions (i.e.) 

Pi = 0.56524, Zi = -0.25, (12) 

which are well inside the ergosphere (see figure 2). 

5 Jet far from the core of M87 

Observations concerning the geometry of the jet far from the core of the 
galaxy M87 have been performed first by Junor et al. (Junor et al. 1999) and 
being interpreted, like with the observations of Doeleman et al. (Doeleman 
et al. 2012), by assuming the validity of the model of Blandford and Znajek 
(Blandford & Znajek 1977). More recent observations, Asada and Nakamura 
(Asada & Nakamura 2012), confirm and specify the data obtained by Junor 
et al Junor et al. 1999). 

The fact that we can interpret the observations obtained by Doeleman 
et al. (Doeleman et al. 2012) in the light of model (Pacheco et al. 2012), 
induces us to ask if it can say something about its geometrical structure far 
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Figure 4: The central curve is the plot of the intrinsic (transverse) radius Pjetint 
(here noted pj) as function of a, deduced from the apparent (radial) radius 
fjetApp = 5.5 observed at the jet base by (Hada et al. 2011). The two in- 
trinsic transverse radii corresponding to the two limits of their margin of error, 
fjetApp ^ [5.1,5.9], are also plotted. For the value a = 0.877004, corresponding to 
the largest ergoregion (see figure 1), the radius at the base of the jet is pj = 3.16372. 



from the core. To answer this question we confront the fittings suggested by 
Asada & Nakamura (Asada & Nakamura 2012) to the geodesies studied in 
the model (Pacheco et al. 2012). 

Asada & Nakamura (Asada & Nakamura 2012) observe an important 
changing of the jet slope after a certain altitude of the order z ~ 10 5 , near of 
the Bondi radius, indicating a divergent tendency. From these observations 
they established the following best fitting laws. For z G [10 2 , 5 x 10 5 ], 



-0.2(^)'orp=0.8(_l) : 



1/6 



and for z e [10 5 ,2 x 10 



71 
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where b\ = 1.73 is called the parabolic type, and b 2 = 0.96 the conical type. 
We plot these curves in figure 5 (see the two straight lines framed by the two 
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geodesies). We seek to confront them with geodesies of the model (Gariel et 
al. 2012; Pacheco et al. 2012). 

To keep continuity with the previous situation, near the core, studied 
in section 4, we keep the parameters values E = 10 6 and p% = 3.1637 to 
look for unbound geodesies, not asymptotic to pi, by numerical integration 
of equation (21) given in (Gariel et al. 2010), for different i.e. inside the 
ergosphere. We plot in fig. 5 an "external" and an "internal" unbound 
geodesies, starting from the ergosphere, obtained for the i.e. {pi = 0.7, 
Zi = —0.4} and {pi = 0.85, z-i = —0.15} respectively, which frame the two 
precedent best fits. The initial parts of these geodesies can be also viewed in 
figure 2, where they frame a set of various geodesies (not represented in the 
figure 5, for the sake of clarity) obtained for intermediate initial conditions, 
all inside the ergosphere. The "internal" geodesic starts above, very close to, 
the geodesic asymptotic to p\ defined in ffl2]) . and later diverges. The flow 
of geodesies located between the external one and the internal one diverge 
also, forming an outer shell, distinct of the perfectly collimated geodesic f|T2|) . 
There is evidence that the jets produced by blazars, are constituted of an 
external envelope, called sheath, and of an interior thin beam (Giroletti et 
al. 2004; De Villiers et al. 2005; Xie et al. 2012). These features appear 
naturally in our purely gravitational model. 

We can see on figure 5 that the "conical" fit (13) is in good agreement 
with the "internal" geodesic (and those close to it), the agreement seeming 
even extend to lower altitudes, up to z ~ 10 3 where the fit is however not 
consistent with the observations. In contrast, the "parabolic" fit ffTBl does 
not seem able to be equated with any single geodesic of the flow. The flow 
intersects the fit. That should be interpreted by noting that the successive 
points emitting the observed radiations belong to the outer portion of the jet. 
Thus, there is an erosion of the width of the jet during its upward motion, 
mainly due to its friction with the ISM. This friction is greater outside, where 
the density of the ISM is greater, because the centrifugal force due to the BH 
rotation drags it outside. Without this erosion effect, the external geodesic 
would remain the external part of the jet which would remain parallel to the 
conical fit from z ~ 10 3 (or even before), as it can be seen on figure 5. 
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Figure 5: Plots, in LogLog scales, of two diverging unbound geodesies: an "ex- 
ternal" geodesic, obtained for the i.e. (pi = 0.7, Zj = —0.4), and an "internal" 
geodesic with i.e. (pi = 0.85, Z{ = —0.15). At their beginnings (see also fig- 
ure 2), they are close to above the perfectly collimated one obtained for the i.e. 
(pi = 0.565256, Zi = —0.25). The two straight lines are the plots of the fits 1 and 
2 obtained by (Asada & Nakamura 2012) in their respective domains of observa- 
tions. We can see that the external and internal geodesies frame the fits. Between 
these two geodesies, a continuous flow of geodesies can exist, some of them being 
represented at their beginnings from the ergosphere in figure 2. All of them are 
external to the geodesic closest to the z-axis with asymptote p\ = 3.1637 ("spine"), 
and diverge from it, shaping the external shell of the jet ("sheath"). 
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